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Field of the Invention 



The invention relates to an ultrasonic imaging system, and/or an ultrasonic 
examination apparatus, for acquiring a sequence of medical images of an artery segment 
using a curved array of transducer elements. The invention further relates to an image 
processing method for processing the images produced by said system and particularly for 
10 displaying an ultrasonic image sequence of an artery segment with indications of arterial 
parameters in function of the cardiac cycle. The invention is used in the field of ultrasonic 
imaging, to provide a cardio-vascular non-invasive medical tool for examining patients 
suspected to present anomalies of arteries and notably anomalies of the aorta such as aortic 
aneurysms. 

1 5 Background of the Invention 

An ultrasonic image processing method for calculating dilation curves related to an 
artery segment is already known from the patent US-05,579,771 (Bonnefous, Dec. 3, 1996). 
This document describes a method for characterizing an artery segment by ultrasonic 
imaging, using an array of ultrasonic transducers that produces a sectional frame, which is 

20 formed by image lines of a number of successive parallel excitation lines extending 

perpendicularly to the artery axis. Said array is coupled to a transmitter/receiver circuit, 
which provides high frequency signals to a signal processing system. Said system determines 
the arterial walls radial velocity and displacement amplitude values and further determines an 
arterial dilation curve in function of location and time. Such a curve is constructed by points 

25 representing the arterial dilation value in the arterial radial direction Z, at a given location 
corresponding to an excitation line along the longitudinal X-axis of the artery, in function of 
excitation instants t, during a cardiac cycle. So, FIG.4C of this document shows, superposed, 
the different dilation curves related to all the excitation lines of an ultrasonic signal 
corresponding to the examined artery segment, said lines being at regularly spaced locations 

30 along the X-axis of the artery. 

A problem is that the cited document relates to an image processing method based on 
image acquisition with ultrasound scanning lines that are perpendicular to the artery axis. 
This corresponds to the use of an ultrasound system for acquiring the ultrasound data with a 
linear array of transducer elements. This kind of system is appropriate for studying a shallow 
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artery and a small segment of artery such as the carotid. This kind of system is not 
appropriate for the study of a deep and thick artery such as the aorta and particularly for the 
study of Abdominal Aortic Aneurysms (AAA). For studying the aorta and AAA, a curved 
array of transducer elements is preferably used. When the ultrasound data are acquired with a 
5 curved array, then the method disclosed in the cited document for calculating artery dilations 
cannot be directly used, since the scanning lines are no longer perpendicular to the artery 
axis. 

In order to early diagnosing aneurysms in aorta, the medical field has a need for non- 
invasive means for providing aorta images together with clear quantified indications of the 
10 aortic dilation. 

Summary of the Invention 

In order to address the problem of finding new diagnosis information for the follow 
up of patients suspected to present Abdominal Aortic Aneurysms (AAAs), it is an object of 
15 the invention to propose an ultrasonic imaging system for acquiring a sequence of medical 
images of an artery segment using a curved array of transducer elements. This system has 
processing means and display means to produce a sequence of images of artery walls with 
specific patterns for visually setting in evidence artery wall anomalies. The system of the 
invention is specially designed to constitute a tool for non-invasive study of anomalies of a 
20 deep artery such as the aorta. 

This ultrasonic imaging system is claimed in Claim 1. 

It is also an object of the invention to propose an ultrasound apparatus having such an 
ultrasonic imaging system. It is a further object of the invention to propose an image 
processing method for displaying an image sequence of deep arteries. In particular, it is an 

25 object of the invention to propose an image processing method for evaluation of parameters 
related to the tension and strain of the aneurysm walls of the aorta. The present invention 
proposes a method developed for AAAs that is specifically designed to provide clinicians 
with information on the motion of the aorta artery walls. This image processing method 
offers the advantage that the aorta wall behavior is made clearly visible together with the 

30 parameters that are useful for the clinician in the study of these Abdominal Aortic 
Aneurysms. 



Brief Description of the Figures 
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Specific embodiments of the invention will be described in detail hereinafter with 
reference to the accompanying diagrammatic drawings; therein: 

FIG.l shows a schematic representation of an aorta and Abdominal Aortic Aneurysm 
(AAA); 

5 FIG.2A is a block diagram showing the main steps of the method of the invention; 

FIG.2Bis a block diagram of an examination apparatus with a viewing system having 
processing and display means for carrying out the method of the invention; 

FIG. 3 illustrates the geometry of the scanning and raw color data storing order; 

FIG.4 illustrates the weights applied to evaluate best depth for intersection between 
1 0 color lines and structures; 

FIG.5 shows an ultrasound image with color lines and structures and the intersection 
between color lines and structures; 

FIG.6 illustrates the principle of motion amplitude correction for continuity between 
cardiac cycles; 

15 FIG.7 represents a color line and a projection line with respective angle from 

reference; 

FIG.8 shows a display of the artery walls and motions for each frame of the sequence; 
FIG.9 shows a user interface to summarize the aorta behavior, with annotations to 
provide information on the meaning of the different lines and sections. 

20 

Detailed Description of Embodiments 

Referring to FIG.l A, Abdominal Aortic Aneurysm AAA is defined by a doubling of 
the normal diameter of the infra-renal aorta A. The heart is denoted by H. The AAA 
abnormality is present in 5% of men aged over 65 years. Rupture of the aneurysm, the most 

25 common complication of AAA, is responsible for about 2% of deaths in men in this age 
group and is the tenth leading cause of death in men in Europe. Since most AAAs are 
asymptomatic until rupture occurs, up to 50% of all AAAs repairs are performed as an 
emergency operation. As the operative mortality for ruptured AAA is around 50%, and only a 
small fraction of patients with ruptured AAAs survive to reach hospital, the overall 

30 community mortality for ruptured AAAs is estimated at over 90%. For this reason, there is an 
increasing interest in the clinical and cost effectiveness of mass screening programs for 
AAAs. Acquired abdominal aortic aneurysms classically are characterized anatomically by 
an unparallelism of the aorta edges, resulting in an expanded and beating abdominal mass. 
The physiopathology consists of a loss of vascular contention, including a risk of rupture. 
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Indeed, the aorta fulfills several haemodynamic functions of blood tissue distribution, 
damping of the pulse wave, etc. The most elementary of these functions is containing high 
pressure blood within the arterial lumen. Arterial wall aneurysmal diseases are characterized 
by partial loss of integrity called dilation or total loss of integrity corresponding to a rupture. 
5 Therefore, in order to early diagnosing aneurysms in aorta, the medical field has a need for 
non-invasive means for providing aorta images together with clear quantified indications of 
the aortic distensibility. Besides, it is important to use non-invasive means instead of invasive 
means because invasive means modifies the aorta pressure, hence the actual aorta 
distensibility. 

10 The severity of Abdominal Aortic Aneurysm (AAAs) is generally clinically estimated 

by considering its maximal diameter. Failure occurs when the wall stress exceeds the strength 
of the artery walls. Whereas operative indications for elective AAA repair are generally 
based on aneurysm size greater than 4.5 to 5 cm in diameter, the most frequently used 
medical approach is watchful waiting, whereby aneurysm diameter is periodically re- 

15 measured to detect expansion to a size warranting surgery of the patient. Now, it is also A 

known that AAAs with a diameter less than 5 cm can rupture. Hence, there is a clear need for 
additional diagnosis information related to dilation. 

The present invention proposes an image processing method to provide aorta 
parameters for the evaluation of the tension and strain of the aneurysms walls. This method is 

20 developed for AAAs and is specifically designed to provide clinicians with information on 
the behavior of the aortic artery walls. This method permits of evaluating automatically, or 
with limited user interaction, and at any time in the image sequence, the position of the artery 
walls, in order to estimate the artery dilations and distensibility. 

Referring to the block diagram of FIG.2A, the processing of an image sequence 

25 comprises the following main steps of: 

1) Acquisition of the image sequence 51: 

This Abdominal Aortic Aneurysm Wall Motion (AAAWM) tool first comprises 
means for acquiring of a sequence of ultrasound images of a segment of artery, for instance a 
segment of aorta, using a linear curved array. Said artery segment has a longitudinal axis and 
30 is represented in grayscale images in FIG.5 or FIG.9. As a matter of example, the processed 
sequence of abdominal aortic aneurysms (AAA) has been acquired using a C5-2 probe and a 
Philips HDI5000 scanner. 

FIG.2B shows a diagram of a medical viewing system 150 according to the invention 
for carrying out the steps of the image processing method described hereafter. The system 
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has means 151 for acquiring digital image data of a sequence of images, and is coupled to 
computer means 153 for processing these data according to this image processing method. 

The data processing device 153 is programmed to implement a method of processing 
medical image data according to invention. In particular, the data processing device 153 has 
5 computing means and memory means to perform the steps of the method. A computer 
program product having pre-programmed instructions to carry out the method may also be 
implemented 

The computer means 153 can be used in the intervention room or near the 
intervention room for processing sequence images. Steps of the present method can be 

10 applied on stored medical images, for example for estimating medical parameters. The 

medical viewing system provides the image data by connection 157 to the system 153. The 
system provides processed image data to display means and/or storage means. The display 
means 154 may be a screen. The storage means may be a memory of the system 153. Said 
storage means may be alternately external storage means. This image viewing system 153 

15 may comprise a suitably programmed computer, or a special purpose processor having circuit 
means such as LUTs, Memories, Filters, Logic Operators, that are arranged to perform the 
functions of the method steps according to the invention. The system 153 may also comprise 
a keyboard 155 and a mouse 156. Icones may be provided on the screen to be activated by 
mouse-clicks, or special pushbuttons may be provided on the system, to constitute control 

20 means 1 58 for the user to actuate the processing means of the system at chosen stages of the 
method. 

The means 151 for acquiring digital image data of a sequence of images may be an 
ultrasound examination apparatus coupled to this medical viewing system 150. This medical 
examination apparatus 151 may include a bed on which the patient lies or another element for 
25 localizing the patient relative to the apparatus. The image data produced by the ultrasound 
examination apparatus 151 is fed to the medical viewing system 150. 

The image processing method of the invention, for forming the Abdominal Aortic 
Aneurysm Wall Motion (AAAWM) tool, is described hereafter more precisely. In the 
following, the terms artery wall border and "structure" have the same meaning and represent 
30 a segmented object. This image processing method comprises further steps of: 

2) Raw color information acquisition and projection in display coordinates 52: 

The segment of artery is further color scanned using the curved array of transducer 
elements. Since the transducer array is curved, the method disclosed in the cited prior art 
cannot be directly used. Scanning with the curved array, in an ultrasound color mode, for 
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example with a Tissue Doppler Imaging (TDI) modality, permits of acquiring ultrasound 
color data related to motion of tissue. The ultrasound raw color data provides the ultrasound 
information for processing the artery wall motion. The ultrasound information is composed 
of the lines or beams of the ultrasound color scanning and, for each ultrasound line, the 
5 estimates of velocities of the artery walls along the ultrasound color lines in depth. Color 
geometry of the scanning and raw color data storing order of the color acquisition is 
represented in FIG.3. Referring to FIG.3, in the method developed for the AAAs, the indexes 
(angle index, depth index) of the ultrasound raw color data must be converted into display 
coordinates (X,Y) in order to match the ultrasound color information with the position of the 

10 artery walls. The equations (la, lb) of the transformation, for the conversion of an index 
(angle index, depth index) in the ultrasound color raw domain into the display domain in 
pixels (X,Y), are presented hereafter. In this FIG.3: C is the scanning center; X0, Y0 are the 
position of C in display coordinates of the display domain; Al is the starting angle; A2 is the 
stopping angle; CL is the current color line; Ref is the Reference Angle =0°. In these 

15 equations: 

StopAngle and StartAngle refer to the stop angle and start angle in radians of the 
ultrasound color information or ultrasound color beam; 

NumAngles is the number of beam angles in the color data; 

8a is the interval in radians between two successive color beams at two respective 

20 angles; 

8depth is the interval in pixels between two successive color estimates on a given 
beam at a respective angle; 

a is the current angle in radian of a current beam corresponding to the angle index in 
the color domain; 

25 depth is the depth in pixels corresponding to the depth index in the color domain; 

(X,Y) is the position of ( angle index, depth index) in the display domain, in pixels. 

8a = (StopAngle - StartAngle) / NumAngles 
5depth = (StopDepth - StartDepth) / NumDepths 
a = StartAngle + angle _ Index * 8a 
depth = StartDepth + depth _ index * Sdepth 

X = X0- sin(cc) * depth (la) 

Y = Y0 + cos(cc) * depth (lb) 

30 3) Evaluation of the intersections between structures and color information 53: 
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The structures are the two internal boundaries of the artery walls that are to be 
previously determined, for example using the method described in the cited prior art. The 
structures are determined using the gray scale images. The structures are reported in the raw 
color data frames. 

5 From now on, the color lines that are used in the display are the display color lines 

whose angles are calculated asabove from the ultrasound color lines. The ultrasound color 
estimates at each depth are calculated and reported on the color lines for display as follows. 

The indexes that correspond to intersections of the color lines with the structures are 
determined in the raw color data frames in the ultrasound color domain. For each frame and 

10 each structure, the points of the structure are associated with their distance to the closest color 
estimate. Each pixel of the structure is associated to the closest line angle in the color 
domain. Then, for all the pixels associated to the same color line angle, the depth is estimated 
as follows: the final depth in the color domain is the gravity center of the depths of the 
considered pixels. The weights W are defined as the inverse of the remainder, denoted by R, 

15 between the closest line angle in the color domain and the line angle between the pixel and 
the scanning center C. For small remainders, the weight function is thresholded as shown in 
FIG.4, which represents the weights W applied to evaluate the best depths for intersections 
between color lines and structures. In FIG.4, 5a is the distance in radians between two 
successive color line angles in radian, and (8a/2) is the maximum remainder R. 

20 A result of the estimation of the intersections between the color lines and the 

structures is presented in FIG.5. The extremities of the small fine lines represent the 
intersection between color lines and the artery wall boundaries that are called reference 
structures. When there is no intersection between a color line and the reference structures, or 
no intersection with only one of the reference structures, the color information of the 

25 corresponding line can not be used to evaluate the artery dilations and thus is not considered 
for further processing. 

4) Velocities averaging 54 

For each structure, the velocity associated with each color line is the result of an 
averaging of a few velocities. The number of indexes selected for the averaging depends on 
30 the width in mm of the wall. The usual value for the wall thickness is for example 1mm. In 
order to limit the influence of noisy data when the velocities are estimated at positions too 
close to the artery border, the averaging can be performed on velocities corresponding to 
positions located more inside the walls. An offset variable has been defined to specify the 
amplitude of the displacement towards the inner parts of the walls. It can also be set to 0 if no 



WO 2004/054448 



PCT/EB2003/005928 



8 



offset is required. A global averaging of the velocities for each of the structures provides a 
mean to estimate the cardiac cycles beginnings for the whole sequence. 

5) Wall Motion estimations 55 

For each color line, the velocities of each structure are integrated over time, for 
5 example over a cardiac cycle, since the cardiac cycles beginnings have been previously 
determined, which permits of determining the cardiac cycle duration. This provides motion 
information of a structure along each color line over time. Since the integration constant 
remains unknown, the motion of a structure is not perfectly cyclic and the amplitude of the 
motion at the end of a cardiac cycle can be different from 0, thus showing a shift S. In order 

10 to provide comprehensible information about motion of the artery walls, a choice is made to 
reset to 0 the amplitude of the motion at the end of each cardiac cycle. In order to keep 
continuity of motion, an affine correction of the data is performed in each cardiac cycle. 
FIG.6 illustrates the principle of motion amplitude correction for continuity between cardiac 
cycles. The amplitude of motion before correction is represented by the curve CI . The 

1 5 amplitude of motion after correction is represented by the curve C2. 

6) Doppler Angle Correction 56 

The amplitude of the wall motion is corrected to compensate for the Doppler angle. 
FIG.7 is a view of the color line and the projection line, with respective angle from a 
reference angle denoted by Ref. In FIG.7, the reference angle Ref is represented by a vertical 

20 line. The angle of the color line denoted by CL is represented as a doted line and its value is 
called a. The angle of the estimated motion direction denoted by MD is represented with a 
bold line and its value is called |3. The angles are taken in the trigonometric orientation and 
are signed. Thus, the resulting Doppler angle between the color line CL and the projection 
line is the difference of a and p. 

25 The corrected motion amplitude is calculated using the following equation (2), which 

gives the motion amplitude correction using the Doppler angle, where WM represents the 
corrected motion amplitude of the measured motion amplitude WMtdi- 



30 The dilation estimation is the result of the difference of motion between the two structures 
(artery wall boundaries) for each color line CL. The dilations are calculated in order to 
provide input data for the interface of the application. The distensibility is the ratio between 
the dilation and the diameter of the artery. 




(2) 



7) Dilations estimation 57 



WO 2004/054448 PCT/IB2003/005928 

9 

8) Motion representation 58 

In order to represent the motion in the images, a choice must be made regarding the estimated 
direction of the motion. In this application, the hypothesis is that the motion of the artery 
walls is perpendicular to the artery principal axis. This is illustrated by FIG. 8, which shows 
5 the motion estimation of each structure represented on the normal to the global artery axis, 
externally to the reference structures. 

9) Display in the frames of the sequence 59 

Referring to FIG. 8, the display provided in each frame of the sequence is limited to 
two types of information. The first type is the structure location. The proximal and distal 

10 walls are represented in color in order to be more easily visualized by the user. For instance 
the two wall structures are represented in the same color. Then, the motion of each wall along 
each color line is represented in a second color to be more easily distinguished by the user. 
The reference line for a null motion is the structure itself, called reference structure, and the 
motion amplitudes are represented starting from the reference structure position. The 

15 representation of the lines of the second color, and direction normal to the artery axis, for 
each motion amplitude, allows of understanding the direction of projection that is selected. 
The lines of said second color are interconnected to represent the overall shape of the motion 
between the lines of the second color. FIG.9 shows a display of the artery walls and motions 
for each frame of a sequence. 

20 10) Display on a dedicated interface 60 

After the processing, the results are summarized on a dedicated interface as shown in 
FIG.9. Figure 9 is a user interface to summarize the aorta behavior, annotated in boxes to 
provide information on the meaning of the different lines and selections. 

The upper left part of the interface presents an echo image, referred to by 10, 

25 corresponding to the user selected frame, combined with the segmentation result for the 

proximal and distal walls and the dilation amplitudes of motion for said proximal and distal 
walls. The selected current color line is shown: 

1 1 are the selected limits of color; 

12 represent the segmentations of proximal and distal walls; 
30 13 are the motions of proximal and distal walls; 

CL is the current color line. 
The middle left part displays curves, referred to by 20, of the maximum amplitude 
and minimum amplitude of the dilations for a given cardiac cycle, as a function of the color 
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lines. The selected current color line is associated to the display of the same color line in the 
echo image 10: 

Max is the maximum dilation per line in the current cardiac cycle; 

Min is the minimum dilation per line in the current cardiac cycle. 
In the lower left part of the display, the amplitudes of the dilations, referred to by 30, 
are displayed as a function of time. It allows of comparing the amplitudes of the dilations 
between different color lines: 

31 is the selected time t in the sequence; 

32 is the current cardiac cycle. 

The upper right part presents the dilation amplitudes, referred to by 40, for the cardiac 
cycle selected by the user: 

41 is the mean dilation for current cardiac cycle; 

42 is the dilation of selected color line for current cardiac cycle. 

The user may have an interaction with clicks of the mouse of the ultrasound system in 
order to select the color line CL in display 10 or 20 or 40; or the time t represented by 31 in 
display 30 of FIG. 10. 



